Two-stage Kondo effect in side-coupled quantum dots: Renormalized perturbative 
scaling theory and Numerical Renormalization Group analysis 
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We study numerically and analytically the dynamical (AC) conductance through a two-dot system, 
where only one of the dots is coupled to the leads but it is also side-coupled to the other dot through 
an antiferromagnetic exchange (RKKY) interaction. In this case the RKKY interaction gives rise 
to a "two-stage Kondo effect" where the two spins are screened by two consecutive Kondo effects. 
We formulate a renormalized scaling theory that captures remarkably well the cross-over from the 
strongly conductive correlated regime to the low temperature low conductance state. Our analytical 
formulas agree well with our numerical renormalization group results. The frequency-dependent 
current noise spectrum is also discussed. 

PACS numbers: 75.20.Hr,74.72.-h 



Introduction. The Kondo effect— in semiconductor 
quantum dots has attracted significant theoretical and 
experimental interest in recent yearsi^i^ii^ In a dot that 
is only weakly coupled to leads, charge fluctuations are 
typically suppressed due to Coulomb blockade^. How- 
ever, if the dot has an odd number of electrons then the 
spin of this electron interacts antifcrromagnetically with 
the spin of the conduction electrons in the leads, and at 
low temperatures it is screened through a Kondo effect. 
The formation of this Kondo state typically leads to an 
enhancement of the conductance at low bias voltages. 

Recently, in the quest for designing multiple quantum 
dot systems with tunable spin control, which can be used 
in spintronics and quantum information processing, dou- 
ble quantum dots have become the focus of interest .•SjI 
In these systems, when both dots are tuned to the sin- 
gle spin regime, an effective spin-spin interaction known 
as the Ruderman-Kittel-Kasuya-Yoshida (RKKY}^ in- 
teraction is mediated between the two dots by the con- 
duction electrons.^ This RKKY coupling competes with 
the Kondo effect in these systems. In the case where 
the two dots are coupled to two separate electrodes, an 
antiferromagnetic RKKY interaction leads to a cross- 
over between the Kondo and RKKY regime o^'^'^^d^d'^ : 
a sufficiently strong antiferromagnetic RKKY coupling 
will lock the spins of the two dots into a singlet, and 
thereby suppress the Kondo effect, while the Kondo ef- 
fect persists for weak RKKY interactions. This picture 
is slightly modified under non-equilibrium conditions: 
Then, even for strong RKKY interaction, the Kondo 
effect is partially restored by the finite bias voltage al- 
lowing for triplet excitations.^^ Part of the rich physics 
has been studied previously in the framework of two im- 
purity Kondc^* and Anderson impurity^ models, and 
the singlet-triplet cross-over has also been recently stud- 
ied experimentally by Craig et al.^ who observed a the 



Kondo resonance for weak RKKY couplings and a split- 
ting of the Kondo resonance for large RKKY interactions. 
Significant theoretical and experimental effort has been 
devoted to the Kondo-RKKY transition in these dou- 
ble quantum dot systems, and related multi-orbital sys- 
tems with two Kondo-screening channels have also been 
studied extensively in recent years by different theoretical 
approaches. 

Rather than studying the usual singlet-triplet transi- 
tion itself, in the present paper we shall focus on an even 
simpler but equally interesting arrangement, where only 
one of the dots is coupled to external leads, but the two 
dots are still side-coupled to each other (see Fig.[T]). Un- 
like the usual four lead set-up, where the cross-over be- 
tween an RKKY and Kondo regimes is governed by a 
quantum critical point corresponding to a finite value of 
the RKKY coupling J,— in this side-coupled system 
the side-coupled spin is always screened for any positive 
J, and a Kosterlitz-Thouless-type quantum phase transi- 
tion and a two-stage Kondo effect occurSfi2ii^ We remark 
that this transition is essentially the two-dot analogue of 
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FIG. I: Sketch of the two side-coupled quantum dots. Elec- 
trons tunnel from only one of them to the leads and they 
interact through an effective exchange interaction, J. 
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the singlet-triplet quantum phase transition found in sin- 
gle dot devices, when the two dot levels are coupled to a 
single conduction electron mode in the leads . 

Although many interesting results have been obtained 
recently for this side-coupled system/^' there is still a 
lack of a detailed analytical and numerical understanding 
of the two-stage Kondo screening processes in the vicin- 
ity of the quantum phase transition. The main goal of 
the present work is to have some more detailed theoreti- 
cal control and understanding of the transport properties 
of this transition. We shall reach this goal by combin- 
ing numerical renormalization group^ and renormalized 
perturbative scahng approaches. As we shall see, the 
latter relatively simple analytical framework is able to 
account for the numerical results over a wide range of 
energy scales, and together with Fermi liquid theory, it 
provides a reliable theoretical framework to understand 
the low energy cross-over. 

We mostly focus on the T-matrix, but also present re- 
sults on the linear AC conductance and the equilibrium 
current fluctuations. Both of these quantities can be ob- 
served experimentally and show features characteristic of 
the two-stage Kondo effect. 

The Model. For the numerical calculations we shall de- 
scribe the two side-coupled dots in Fig.[T]by an Anderson- 
like model. The two isolated dots are described by the 
Hamiltonian 



H 



DD 



+ JS1S2, 



(1) 



where i — 1,2 labels the two dots, Ni — Acr'^ia 
is the number of electrons occupying dot i and Si = 
(1/2) dl^CTCTCT'dj^/ is their spin. Each dot is sub- 
ject to a charging energy, Ui ~ U2 = U = Eq, and the 
two dots are coupled by an exchange coupling, which 
is assumed to be antiferromagnetic, J > 0. The rigi 
in Eq. ([1]) denote dimensionless gate voltages that set 
the occupation numbers, {Ni). In the rest of the paper, 
we restrict ourself to the case of particle-hole symmetry, 
rigi = ng2 = 1. However, this assumption is expected to 
be unimportant as long as Ugi w ng2 ~ 1. 

The coupling of dot 1 to the leads is modeled by the 
usual tunneling Hamiltonian, 



(2) 



a=L,R e,(T 



Here, Vl and Vr denote the tunneling amplitudes to the 
left and right leads, respectively, and cj^^g. creates an elec- 
tron in lead a = L, R with spin a and energy e. This 
tunnel coupling leads to a broadening of the level on 
dot 1, the width of which is given by F — F^, -I- Fi^ = 
27r(y^PL + V^gii}, with Ql/r the density of states in the 
leads. 

AC conductance and noise: Our main goal is to de- 
termine AC transport properties in the linear response 
regime and see how the two-stage effect appears in these 
quantities. Fortunately, since there is no charge transfer 



between the two dots, the derivation of Ref. carries 
over to our case, and the real part of the optical conduc- 
tance is simply given by 

^'H^r^E / du;'[f{u;'-u;)-f{cj'+u;)]T,{cu'), (3) 

where T(uj) is "the transmission probability" at energy 
uj, and 



Go — 



(4) 



denotes the maximum conductance through the dot. The 
transmission coefficient Ta-{uj) appearing in Eq. ^ can 
be expressed as 



Tcy{uj) = -F Ini Giiait^) 



(5) 



where we have introduced the retarded Green's function 
on dot 1: Giiait) = -i6'(i)({di^(t), d^}). This Green 
function can be computed accurately using numerical 
renormalization group methods, but as we shall see, 
substantial analytical progress can be made by a special 
version of renormalized perturbation theory. Although 
Eq. ([5]) is valid at any temperature, in the following we 
shall focus our attention to the high frequency regime, 
uj ^ T, and set the temperature to zero, T = 0. 

The noise spectrum of the device is also of experimen- 
tal relevance. This is defined as 



Cico) 



dte'"*[(/(0)/(i)) - 



(6) 



At equilibrium, this is simply related to the linear con- 
ductance G'{uj) by the fluctuation-dissipation theoremj^^ 



2huj 



exp{huj/kT) — 1 



(7) 



This formula simplifies further at T = temperature 
to C{uj) = 2h\uj\ G'{ijj) 0{—Lo). Clearly, to determine 
the equilibrium noise spectrum of the device and its AC 
conductance only the Green's function of the first dot 
needs to be determined. 

Renormalized perturbative scaling: Before presenting 
our numerical results, let us reach some analytical un- 
derstanding of the physics of the side-coupled dot in the 
limit J ^ 0. In this regime, "two-stage Kondo screening" 
takes placeii: In the first stage, the spin of the dot 1 gets 
screened below the Kondo temperature Tk ~ D e~^^^^ , 
where the high energy cut-off D denotes the effective half 
bandwidth of the conduction electronsi^ Clearly, for the 
Kondo effect to take place J ^ Tk is required, otherwise 
the two spins are locked together to a singlet before the 
Kondo effect can take place. Then below the Kondo scale 
Tk the electron on the first dot is dissolved in the con- 
duction electron sea of the leads, and presents an effective 
fermionic bath for the electron on the second dot. Since 
the coupling between the two dots is antiferromagnetic. 
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FIG. 2: Top: leading logarithmic contribution to the self en- 
ergy of the first dot's Green's function. Solid lines represents 
the propagator G^i of quantum dot 1, given by Eq. ((8)1, while 
dashed lines denote the pseudofermion propagator associated 
with the side-coupled dot's spin, 82- Full squares stand for 
the full leading logarithmic vertex function, given in the lower 
part of the figure. Crossings of dashed and continuous lines 
correspond to an interaction through the renormalized RKKY 
interaction, Jeff — J/z, given by Eq. Q. 



another Kondo effect shall take place at a much smaller 
energy scale, T*, where the spin of the second dot is also 
screened. Our aim is to understand the formation of this 
second Kondo singlet in detail. 

We first observe that in the regime of interest, ui <C Tk 
the Matsubara Green's function of dot 1 can be approx- 
imated by the resonant level expression,— 



iun + i TKsgn{uJn) 



(8) 



where z — 



denotes the quasiparticle weight at the 



Fermi energy, and Tk = zT — c Tk is an energy of the 
order of the Kondo temperature, Tk- The precise value 
of the universal constant c relating Tk and Tk depends 
on the definition of Tk- Throughout this paper we shall 
define Tk as the half-width of the transmission T(a;). 
Then from fitting the NRG data we get c « 0.5. Note 
that the Lorentzian representing the Kondo resonance 
has a very small spectral weight, 2; ^ 1, and most of the 
spectral weight goes to the Hubbard peaks. 

Fermi liquid theory and the basic principles of renor- 
malization group also imply that the exchange interac- 
tion between the dot spins is renormalized by the same 
z-factor, i.e., in the regime oj -C Tk the effective RKKY 
interaction reads 



-ffRKKY 



J 



Si S2 



(9) 



In a first approximation one would think that, J = J, 
however, the fact that ImGdi has a large logarithmic tail 
above Tk leads to a slight renormalization of this rela- 
tion. From a fitting of the numerical data shown later, 
we obtain the approximate relation J Ri 1.1 J. 

We are now in the position to develop a perturbative 
scaling theory in the small coupling, J. To do this, we 



used Abrikosov's pseudofermion representation to com- 
pute the second order self energy and vertex corrections, 
shown in in FigO The dimensionless vertex function is 
given by the following expression: 



-f{uj) = g{uj)T{uj) = giuj)J + {g{u)jflog 



Tk 



(10) 

where Tk — c Tk is the effective width of the Kondo 
resonance, u is the energy of the incoming electron, and 
g{uj) denotes the rescaled effective density of states of dot 
1, serving as a many-body reservoir for dot 2, 



g{uj) = giuj)/z 



T 



K 



(11) 



Eq. (llOp is only of logarithmic accuracy, irrelevant terms 
of order oj/Tk have been neglected. 

The Kondo temperature Tk ~ Tk appears in Eq. (fTO|) 
as a high energy cut-off. We can therefore perform a scal- 
ing transformation by reducing this cut-off, Tk — > T'j^, 
and requiring the invariance of the vertex function for 
frequencies uj < Tk at the same time. This transforma- 
tion sums up all leading logarithmic diagrams and leads 
to the following scaling equation 



d{g{Lo)J) 
dl 



(12) 



with the scaling variable defined as I = \og{TK /T'j^). In- 
tegrating this differential equation up to Z = \og{TK /oj), 
one obtains the dimensionless vertex function in the lead- 
ing logarithmic approximation: 



i{^,Tk) = 



g{uj)J 



1 + g{uj)J logi-uj/TK) 
This equation can be rewritten with a little algebra as 

1 



(13) 



i{.^,Tk) 



1^ log 4^ log 



T* 



with the second scale T* defined as 



T* = Tk exp(-^ Tk/J) 



(14) 



(15) 



This scale also appears in the slave boson approach of 
Ref. [l3, however, the latter approach does not account 
for the logarithmic corrections, which are of our main 
interest here. Clearly, the dimensionless vertex diverges 
at an energy, w « T*, implying that the effective RKKY 
interaction becomes dominant below this scale and di- 
verges in the limit w ^ Oi^ In other words, below T* a 
singlet is formed from the two dot spins, and the scale T* 
can thus be viewed as the effective singlet-triplet binding 
energy. Note that the spin-singlet is always the ground 
state for any arbitrary RKKY coupling J > 0. 

The point J = is special: It separates the ferromag- 
netic phase (J < 0) from the antiferromagnetic phase 
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discussed so far (J > 0). While in the antiferromagnetic 
phase the two dot spins are locked into a singlet, in the 
ferromagnetic phase the spins of the dots are bound to a 
triplet, which is then partially screened by the lead elec- 
trons, and correspondingly, the ground state is a doublet 
for J < 0. Thus for J < the side-coupled system has a 
residual entropy. At the critical point, J = 0, the scale 
T* diverges exponentially, corresponding to a Kosterlitz- 
Thouless phase transition between these two statesJ^ii^ 
The second order self-energy correction to the retarded 
Green's function Gj^^ simply gives the expression 



where S=l/2. Note that this correction scales as ~ 
J"^ /{Tk z) and, although it looks to be very large at 
a first sight, it is actually small compared to (G^i)"^ ~ 
Tk/z as long as J <^ Tk- In the leading approximation, 
this self-energy results in the following Green's function: 



K 



i Tk 



S(S+1) 



(17) 



J' 



EH = 5(5 + 1) — 



1 



4z u + iTi 



(16) 



K 



Summing up the leading logarithmic corrections to the self-energy simply amounts to replacing J in Eq. 
"/(lu)/ g(uj), and thereby results in the following transmission coefficient. 



m by 



Tcriio) = -Im 



(fl - \S[S + iWiLo^ + nh^co)) + ^fK (fl + i5(5 + iWiLo^ + nhHc)) 



(18) 



The logarithmic corrections hidden in 7 result in the for- 
mation of a dip in the spectral density of dot 1, corre- 
sponding to a suppression of the transmission coefficients 
at energies uj ^ T* Tk- This dip in To-(w) also implies 
the appearance of a dip in the AC conductance discussed 
later, and is a clear signature of the formation of a sin- 
glet ground state. Physically, it is a consequence of the 
fact that electrons promoted from one side of the device 
to the other must first break up the singlet of energy T* 
formed by the two dot spins. Clearly, electrons of energy 
Lj <T* are not energetic enough to break up this singlet 
and therefore their transport is suppressed. 

The logarithmic approximation breaks down below 
uj T* . There a Fermi liquid is formed and T^{uj) scales 



Taiio) 



■(r*)2 



[u « T*) 



(19) 



The constant a vanishes for electron-hole symmetry and 
remains typically small unless electron-hole symmetry is 
dramatically broken, while the coefficient 6 is a number 
of the order of unity. 

Comparison with NRG. The transport properties of the 
side-coupled quantum dot system can also be studied nu- 
merically by Numerical Renormalization Group (NRG) 
methods. The transmission coefficient T(lj) through dot 
1 obtained from NRG at different RKKY couplings J is 
plotted in FiglJl In all figures we compensated for a 6% 
loss of the spectral weight. The large resonance is a man- 
ifestation of the Kondo effect displayed by the first dot, 
and the sharp dips in the transmission at w « are due 
to the formation of the singlet state below the energy T* . 

The numerically obtained transmission coefficients are 
compared to the analytical formula in Fig. 01 The 



perturbative expression agrees well with the numerical 
results over a wide range of energy scales between T* 
and Tk- It is interesting to observe the deviations above 
Tk, where the simple Lorentzian approximation we made 
for G^^^ fails to account for the fact that the resonance 
on dot 1 is also a Kondo resonance. This Lorentzian 
approximation thus completely neglects the large loga- 
rithmic tails for lo > Tk- In the inset we also show the 
scale T* as extracted from our fits as a function of 1/J. 
The extracted scales compare very well with the analyt- 
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FIG. 3: Transmission coefficient through dot 1 at zero tem- 
perature for different RKKY couplings, J. The energy unit 
is the half bandwidth, D = 1. The Anderson model's pa- 
rameters were U = 1, e^i = £^2 ~ —0.5 , Fil = Vir — 0.1, 
resulting in a Kondo temperature of Tk ~ 0.0055. We used 
a discretization parameter A = 2. 
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ical expression, Eq. ([151) . indicated by the solid line. 




FIG. 4: Fit of the numerically obtained transmission coefh- 
cient T{uj) by the perturbative expression, (|18|) for various 
values of J. For all fits we used Tk = 0.5Tk = 0.003. The 
inset shows the scale T* extracted from the fit, as a function 
of 1/J. 

The two-stage Kondo effect and the K-T transition can 
also be observed in the AC conductance of the two-dot 
device or the noise spectrum. Such measurements have 
been indeed done in recent experiments of high-frequency 
current fluctuations^, though the linear AC conduc- 
tance measurements in the relevant frequency regime 
are very difficult due to background currents from the 
capacitors According to Eq. ([3]), the real part of the 
conductance through the device can be computed from 
the transmission coefficient through a simple integration. 
The resulting curves are displayed in Fig. O 

It is interesting to remark that the low energy cross- 
over of the transmission coefficient and the conductance 
is described by universal cross-over functions for lu, T* <C 
Tk- The conductance, e.g., is approximately given by 



G'iu;) « Gogiuj/T*) 



(20) 



in this regime. Here the scaling function g depends 
somewhat on electron-hole symmetry breaking, but in 
case of electron-hole symmetry it is completely univer- 
sal. Then for very small frequencies it scales to zero as 
g{u)) « 0.06 {uj/T*)^, while at high energies it approaches 
1 logarithmically, g{uj/T*) « 1 - a/ log^(w/T*). This 
universal cross-over function can be extracted from the 
NRG results, and is displayed in Fig. [B] The transmis- 
sion coefficient T^iuj) displays similar universal scaling 
properties. 

Finally, let us discuss the current noise at zero tem- 
perature, plotted in Fig[71 At T = temperature C[ijj) 
has only weight for a; < 0. For T* < ~uj < Tk we 
have C(a;) cc \io\ corresponding to the Fermi-liquid prop- 
erty of a perfectly transmitting quantum dot^SS. However, 
for — w < T* this behavior crosses over to a power law 
scaling, where in case of electron-hole symmetry one has 




FIG. 5: AC conductance normalized by Go for the same 
parameters as in Fig. |31 



C{uj) oc |a;p/(T*)^. This behavior is somewhat modi- 
fied once electron-hole symmetry is broken. Then the 
conductance remains finite even in the w — > limit, and 
correspondingly, another cross-over may take place from 
the C{uj) cx |wp regime to a linear regime, C(w) oc |a;| at 
some energy T** ^ T*. 

Conclusions. We provided a complete analytical and 
numerical analysis of the AC transport properties and 
the two-stage Kondo screening in side-coupled quantum 
dots. Our analytical results were based on renormalized 
perturbation and scaling theory, and they agree well with 
the NRG result over a wide frequency range. We also de- 
termined the linear AC conductance and the equilibrium 
current noise which can both be measured experimen- 
tally and reflect the two-stage Kondo effect and the K-T 
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FIG. 6: G'{u) as a function of u>/T* for various values of J. 
Below the Kondo scale all curves collapse to a universal curve. 
Inset: small frequency part of this universal curve on a linear 
scale. We also show the quadratic behavior characteristic of 
electron-hole symmetry, determined from a log-log fit of the 
curves. 
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transition. We also computed the universal cross-over 
functions that describe the emergence of the triplet state 
at the energy T* . 
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FIG. 7; Equilibrium current noise C{uj) in case of electron- 
hole symmetry for uj < 0. C{ui) exhibits linear dependence 
on u} for T* < < Tk, that crosses over to a cubic behavior 
below T*. Note that C{ui) — for u > 0. The parameters 
are the same as in FigO 
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